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The  models 

Having  performed  a  study  on  the  orientation  density  function  (o.d.f.)  of  small 
spheroidal  aerosol  particles  in  a  general  field  of  an  arbitrary  strength,  (1,2)  it 
became  possible  to  treat  the  orientation  problem  in  a  turbulent  medium.  To  this  end, 

A — ' 

two  interconnected  physical  models  were  applied.  In  the  first.  The  Realizations 
Model,'H*~ it  was  assumed  that  the  turbulent  particle  field  constituted  an  ensemble  of 
an  infinite  number  of  realizations,  j,  each  one  of  which  is  characterized  by  one  set 
of  the  o.d.f.,  fU)  (x,t),  values.  The  latter  was  taken  in  that  model  to  essentially 
coincide  with  the  previously  found  solution  of  the  Fokker-Planck  equation  (-l-j^'f'in  the 
field  of  the  realization,  r 


3F(^3t  - 
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where  o  is  the  rotational  velocity  of  the  particles  and  c*.  is  their  rotational 

r  '  C  '  * 

Peclet  number  defined  by  W,/  De  ,  W0  being  a  typical  component  of  the  (fluid) 
gradient  tensor,  i>e  an  effective  rotational  diffusion  coefficient  and  R  the 
particle  aspect  ratio;  is  a  gradient  component  d  u-t  /"dXj,  (  - 

a  component  of  the  fluid  velocity).  The  use  of  the  solution  is  based  on  the  estimate 
that,  even  for  the  highest  (Kolmogoroff)  frequency  component  of  u-t  ,  the  rotational 
Reynolds  number  for  the  studied  particles  is  small  enough  to  render  their  motion 
(quasi)  stationary. 
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In  the  secqnd,  "irhe  Micro- turbulence  Model /"a  relation  between  the  so-called 
turbulent  rotational  diffusion  coefficient  of  the  particles  and  the  physical 

characteristics  of  the  fluid  field  was  supplied.  This  coefficient  and  the  Brownian 
rotational  diffusion  one,  ,  compose  the  effective  diffusion  coefficient 


3>e  =  I>t 
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Expressing  the  stochastic  quantities  of  the  system  as  F  =  F  +  F',  col  =  o.  +  o'.  f 
W=  and  u.t=  u;,  +  U'L  ,  it  could  be  shown  that  the  realizations' average 

of  F,  w^>  it, ,  JO,  ,  R  ),for  a  space-time  point  is  actually  F  (wlK,r>ej«?) 

for  that  point  and  that 
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(3) 


The  particles  considered  were  taken  to  be  much  smaller  than  the  Kolmogoroff  scale; 
so,  the  turbulent  rotational  diffusion  coefficient  itself,  rr>  t  ,  was  assumed 
in  the  second  model  to  depend  on  the  randomizing  action  of  the  turbulent  pressure 
fluctuations  at  the  particles'  surface  which  arise  from  the  (Kolmogoroff)  micro¬ 
turbulence. 

From  dimensional  analysis  it  was  obtained  that 
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where  €  ,  the  turbulent  dissipation  energy,  is  given  for  an  homogeneous  steady  field 
by  (3) 
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As  the  turbulent  gradient  tensor  W';.*  is  not  completely  known  in  many  fluid  systems, 
it  was  necessary  in  the  present  study  either  to  supplement  the  missing  data  by  the 
results  of  numerical  simulations  or  (reasonable)  guesses  or  both.  In  cases  of  interest, 

i 

there  was  conducted  a  parametric  investigation  in  which  the  effect  of  structure  changes  'v- 
in  the  gradient  tensor  on  the  o.d.f.  was  tested. 

The  numerical  simulations  of  the  turbulent  field  were  carried  out  according  to  the 

i 

method  of  Wang  and  Frost  (4)  where  u\  =  u-  (t)  is  found.  In  this  method,  however,  ; 
only  the  three  components  W'u  (i  =  1,2,3)  could  be  acquired  due  to  the  applicability  f1 
of  the  relationship  known  as  the  Taylor  hypothesis,  viz. 


a  u.  /  -d  t  = 


a  cl,  /  ax 


to  those  components  alone  . 

r 

The  rest  of  the  gradient  components  were  extracted  either  from  the  eyperimental  findings  [ 
of  Klebanoff  (5)  for  the  studied  situations  of  a  turbulent  flow  over  a  flat  surface  as 
in  the  atmospheric  boundary  layer  or  from  those  of  Wignansky  and  Fiedler  (6)  for  the  • 
studied  situation  of  a  free  round  turbulent  jet,  or  from  parametric  checks  as  mentioned  |j 
above.  The  experimentally  acquired  gradient  components  are  given  in  the  form  of 
a  w'i,,  =  (  )  1/2  and  so  were  taken  the  rest  of  the  components.  The 

(turbulent)  fluid  was  assumed  to  be  incompressible,  viz.  W  a  =  0,  and  it  was  set 


a  =  a  w'„ 
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or  with  reversed  signs. 

Obviously,  since  in  addition  to  the  average  value  of  F  ,F,  or  its  maximal  value  F-^  , 
some  measure  of  the  spread  of  the  o.d.f.  values  is  desired,  a  procedure  was  adopted  in 
which  this  (F)  function  was  calculated  for  the  following  three  realization  fields: 


Finally,  to  account  for  any  uncertainties  in  the  turbulent  field  data,  the 
(determining)  effective  Peclet  number  was  parameterized  by  a  numerical  factor 
to  be 

<*«.  ot0  =  W0  /  “d  e  (9) 

where  W.  wa  W0’  ,  W0‘ [=  (a  ]  being  a  typical  turbulent  component. 
Results 

Two  physical  situations  were  studied: 

1.  The  near-ground  atmospheric  boundary  layer,  in  which  only 

4 

U,  (2)  =  ( u */  ko )  6r»  ( I  -+  Z  IZ0  )  (10) 

and 

Wl3  =Du,l'6x  =.  u*/  K,(Z-Z0)  (11) 

essentially  exist,  z  being  the  elevation  above  ground,  (Z,  the  mean  horizontal 
velocity,  u*  the  friction  velocity,  von  Karman  Constant  and  z  the  roughness 


2.  The  turbulent  round  free  jet,  for  which  the  Schlichting's  solution  for  the 
average  velocity  components  (7)  was  used.  In  this  solution,  the  (molecular) 
kinematic  viscosity  vis  replaced  by  a  virtual  (turbulent)  kinematic  viscosity 
given  in  the  equation 


O.  16  1 
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where  J'  (=J/$  )  =  1.59  bl/a  uc  ,  J'  is  the  kinematic  momentum  of  the 
jet.  bJs  is  its  half-width  and  u0  is  the  average  fluid  velocity  along  its  axis 
(b?i  =  cx  j  in  which  c  is  0.63  to  0.79,  as  experimentally  found). 

In  the  boundary  layer  situation,  both  the  case  of  a  "weak  turbulence",  where 
w„'  We  ,  and  a  "strong  tubulence",  where  W'0  »  w„  ,  were 

investigated. 

For  the  first  case,  numerical  simulations  were  employed  while  for  the  second  one 
Klebanoff's  experimental  data  (5)  supplemented  by  simulations  was  applied.  Also,  in 
the  latter  case,  the  version 


or 
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(13a) 
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was  taken.  The  calculations  in  both  cases  were  carried  out  with^=0.15  cm2/sec., 
h  (thickness  of  the  atmospheric  boundary  layer)  =  103m,  z0  =  0.1m,  k0  =  0.4,  r°B  = 
1  sec-*,  and  e  deduced  from  Eq.  (5).  The  calculation  time  for  F  was  greater  than 
the  relaxation  time  of  the  particles,  and  A0-A<V  =  Tr/i2  or  TT  /  24 


In  the  turbulent  jet  situation,  the  ratio  between  the  various  values  of 
was  obtained  (6)  through 


where  C2=  1  +  exp  (-200r2),  r  being  the  radial  distance  within  the  jet. 

For  both  situations,  the  absolute  values  of  a  WlK  were  deduced  by 

normalization  according  to  Eq.  (5).  Values  of  F,  and  Fm,  together  with  the 
deviations  of  F  for  fields  (2)  &  (3)  (above),  were  calculated  as  a  function  of  height 
above  ground  (and  hence  e  )  in  the  boundary  layer  situation,  the  aspect  ratio  of  the 
particles,  R  (R  >  1  for  fibers  and  R<  1  for  platelets),  and  the  parameterization 
factor  ol q  .  However,  only  the  following  typical  figures  are  presented  here: 


Figure  1.  The  average  o.d.f.,  F  ,  vs.  0  for  the  atmospheric  boundary  layer, 
"weak  turbulence",  with: 

z  =  2.5  to  20m  (1)  2.5m,  (2)  5m,  (3)  10m,  (4)  20m, 

R  *  10,  -  tt  /2, 

»e  =  1(  wi3(=W0  )  according  to  Eq.  (11),  V,.1  =  a  V  \3  .  • 

Solid  lines  are  simulation  results  (254  realizations,  usually);  points  are  values 
calculated  from  the  field. 


Figure  2.  The  maximal  o.d.f.,  Fm,  vs.  height  for  the  atmospheric  boundary  layer, 
"strong  turbulence".  R  =  50;  e  decreases  with  height  according  to  Ball  (8). 

(1)  wj?  ,  (2)  wfj  ,  (3)  vf?  ;  dashed  lines  are  cases  where  terms  of  Eq.  (7) 

are  taken  with  opposite  signs;  points  relate  to  combination  of  Eq.  (13b). 

(Line  1  essentially  coincides  with  the  random  distribution  one  ) 


Figure  3.  The  maximal  o.d.f. ,  Fm,  vs.  particle  aspect  ratio  R  for  the  atmospheric 
boundary  layer,  "strong  turbulence". 

z  =  5m,  €  =  230  cm^/sec-^ 

(1)  ¥<}),  (2)  W(^,  (3)  W(^ 

(Line  1  essentially  coincides  with  the  random  distribution  one  ) 


Figure  4.  The  maximal  o.d.f.  Fm  vs.  the  parameter  o„  f0r  the  atmospheric 
boundary  layer;  "strong  turbulence". 


R  =  50 


i  R  =  -.02  -- 


-  ;  z  =  2.5m;  e  =  5  cm2/sec3. 


Figure  5.  The  maximal  o.d.f.,  Fm,  vs.  particle  aspect  ratio  R  for  a  free  turbulent  jet 
Location:  Xj=  2,  x^=  x3=  0.05  (see  ref.  1), 

"0  =5  sec-3,  £  =  250  cm^/sec3 

(1),(1,):-W(^;(2),(2I)  (3)  ,(3* )  :W(3j*  -  for  (l‘),  (2'),  (3* )  terms  of  Eq.  (7) 

are  taken  with  opposite  signs. 


(Line  1  essentially  coincides  with  the  random  distribution  one  ) 


Figure  6.  The  maximal  o.d.f.,  Fm  vs.  the  parameter  &<,  for  a  free  turbulent  jet 
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1.  The  average  o.d.f.  in  a  weak  turbulent  field  shows  structured  (preferred) 
orientation. 

2.  The  maximal  (and  average)  o.d.f.  in  a  strong,  commonly  occurring,  turbulent 

field  of  the  average  realization  essentially  coincide  with  the  random 

distribution  while  the  spread  of  the  values  of  the  function  between  fields 
(21  (31 

W'^  and  Wv^  is  quite  significant.  This  spread  of  values  may  have  practical 
connotations. 


3.  The  values  of  Fm  and  its  deviations  increase  with  the  parameter 
Free  turbulent  jet: 

Conclusions  2  and  3  of  the  former  situation  apply  here  too. 


as  expected. 
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